VC>2' a two-fluid incoherent metal? 
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We present ab initio LDA+DMFT results for the many-particle density of states of VO2 on the 
metallic side of the strongly first-order (T-driven) insulator-metal transition. In strong contrast to 
LDA predictions, there is no remnant of even correlated Fermi liquid behavior in the correlated metal. 
Excellent quantitative agreement with published photoemission and X-ray absorption experiments 
is found in the metallic phase. We argue that the absence of FL-quasiparticles provides a natural 
explanation for the bad-metallic transport for T > 340 K. Based on this agreement, we propose 
that the I-M transition in VO2 is an orbital-selective Mott transition, and point out the relevance 
of orbital resolved one-electron and optical spectroscopy to resolve this outstanding issue. 

PACS numbers: 71.28+d,71.30+h,72.10-d 

Metal-insulator transitions (MIT) in early-transition 
metal oxides (TMO) are classic examples of correlation- 
driven Mott-Hubbard transitions [1]. Strong multi- 
orbital correlations in the real crystal structure (RCS) 
have turned out to be indispensable for a quantitative 
description of the MIT in these cases [2] . 

A detailed understanding of the above cases is still 
somewhat elusive. In particular, it is still unclear 
whether the MIT is orbital selective, i.e., whether differ- 
ent orbital-resolved densities-of-states (DOS) are gapped 
(driving the Mott insulating state) at different values 
of U, U' (defined below) or whether there is a sin- 
gle (common) Mott transition at a critical interaction 
strength [3,4]. Conflicting results, even for the same sys- 
tem, and within the same (d = 00) approximation [4] 
have been obtained, necessitating more work to resolve 
this issue. 

In this communication, we address precisely this issue 
in vanadium dioxide (VO2). Using the state-of-the-art 
LDA+DMFT [5] technique, we first show how excellent 
quantitative agreement with the full, local one-electron 
spectral function is achieved using LDA+DMFT(IPT). 
We then build on this agreement, claiming that the MIT 
in VO2 is orbital selective, and represents a concrete, 
ab initio realization of the two-fluid scenario for MIT in 
strongly correlated systems. 

V0 2 shows a spectacular MIT at T = 340 K from a 
low-T monoclinic, Mott insulating phase with spin dimer- 
ization along the crystallographic c-axis to a local mo- 
ment paramagnetic metallic phase. This T-driven MIT 
has attracted intense scrutiny: because of the additional 
complication of the antiferroelectric displacement of the 
VOq octahedra in a correlated system, an unambiguous 
characterization of the MIT (relative importance of Mott- 
Hubbard versus Peierls dimerization) is somewhat diffi- 
cult. Various observations have been cited in support of 
both scenarios [1], and theoretical models detailing the 
importance of these effects have been proposed [6,7]. 

One-electron spectroscopies constitute a reliable fin- 
gerprint of the changes in the single-particle spectral 




FIG. 1. Photoemission spectrum of VO2 across the 
metal-insulator transition, from [11]. Top: raw data; bottom: 
with the 0-2p background subtracted in a way indicated by 
the region below the solid (300 K) and dashed (350 K) lines. 

function (and hence those in the correlated electronic 
structures) across a phase transition. Many papers have 
studied the MIT in VO2 by photoelcctron- and X-ray ab- 
sorption spectroscopy (PES and XAS) [8-10]. A perusal 
of the available literature requires care to pinpoint robust 
features (weakly dependent on incident photon energy); 
in what follows, we use only those published data which 
satisfy this criterion in order to make comparison with 
theory. Usually, the experiments have been carried out 
with low incident photon energies (50 — 70 eV) except in 
Ref. [11] (1200 eV). However, the main features, namely: 

(i) a broad step-like feature crossing Ep in the R phase, 

(ii) the "lower Hubbard band" centered at —1 eV, and 

(iii) the huge O 2p-band contribution from —2 eV, seem 
to be roughly consistent in all these measurements. Upon 
subtraction of the 0-2p contribution, the correspond- 
ing "Vanadium" t2 g DOS shows noticeably more spectral 
weight at lower energies relative to that at —1 eV (Fig. 1). 
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In this context, great care has to be taken, as, in view of 
the dynamical transfer of spectral weight over large en- 
ergy scales characteristic of a correlated system, it is not 
a priori clear whether the subtracted "background" con- 
tains significant intrinsic strong correlation (continuum) 
contributions. Additional subtraction due to surface con- 
tributions to the incoherent part [11] results in noticeably 
more spectral weight at low energy compared to that at 
-1.5 eV. 

With these caveats, a consistent picture of the changes 
in the actual electronic structure across the MIT seems 
to emerge. First, the transition is clearly driven by the 
dynamical transfer of spectral weight from high to low en- 
ergy over a large energy scale of 0(2 eV), characteristic 
of a Mott-Hubbard scenario. Second, in the metallic (R) 
phase, there is no semblance of a Fermi liquid-like peak at 
low energy: rather, the low-energy feature is more-or-less 
incoherent. This is in good accord with the high resis- 
tivity in the (metallic) R phase, as also with the linear- 
in-T dependence without saturation upto 900 K [12]. It 
is natural to postulate a common electronic scattering 
mechanism leading to such observations. In the insulat- 
ing (I) phase, clear evidence of insulating gap formation 
is visible in PES. From Fig. 3 of Ref. [11], we estimate 
the gap value Eg — 0.2 — 0.3 eV, consistent with that 
extracted from optical data [13]. 

The XAS data probe the unoccupied part of the full 
spectral function, and show a number of interesting fea- 
tures in VC>2- First, we notice that the XAS intensity is 
much larger (for < oj < 2.0 eV) than the XPS intensity 
(for —2.0 < u> < eV). Second, the broad low-energy 
"step-like" feature seems actually to be the lower part of 
the intense feature in XAS pulled below Ep, and is not 
a thermally smeared correlated Fermi liquid resonance. 
As mentioned above, this fits in nicely with the trans- 
port data; in any case, we do not expect a quasicoherent 
FL scale to survive beyond very low temperatures in a 
correlated system close to a MI transition. 

Beginning with the LDA bandstructure of VO2 in the 
monoclinic crystal structure corresponding to the insulat- 
ing phase [14], with subtraction of terms to avoid double- 
counting of interactions treated on average by LDA, the 
one-electron Hamiltonian is, 

H LDA = £ a/3( k ) C La C k/3a + ^ e iL?W, (1) 

ka/3(T iota 

where e° Qcr = e iaa ~U{n aS - ±) + \oJ H (n aa - 1). Here, 
U, U', Jh are the Coulomb interactions in the ti g shell, 
and are defined below. The lowest lying xy (in the 
M notation) band is the most heavily populated, while 
the yz ± zx bands are less populated. The bonding- 
antibonding splitting [14] is clearly visible in the re- 
sults, which also show that the total bandwidth is about 
W ~ 2.3 eV, contrary to much smaller previous [7] esti- 
mates based on model calculations. Sophisticated LDA 



calculations do yield a miniscule, almost zero-gap, Peierls 
insulator [6], but obviously do not reproduce the Mott 
insulator. As emphasized by Zylbersteyn et al. [7], it is 
inconceivable that the antiferroelectric distortion alone 
could open up a charge gap {Eq c± 0.6 eV), requiring 
proper incorporation of correlation effects. 

Including multi-orbital correlations, the full Hamilto- 
nian is 

ia iaftaa' 

- Jh ^2 ' ^P- (2) 

This multi-orbital Hubbard model is solved using 
LDA+DMFT(IPT) [3]. Accurate low- and high energy 
behavior is ensured by requiring strict numerical com- 
pliance with the Friedel-Luttinger sum rule (low-energy) 
and the first three moments of the DOS (high energy). 
In our approach, the MI transition is driven by the large 
change in the dynamical spectral weight transfer accom- 
panying small changes in the renormalized antiferroelec- 
tric displacement in the correlated picture. Hence the 
transition is of the Mott-Hubbard type, but the antifer- 
roelectric coupling also plays an important role: in fact, 
it acts like an "external field" in the orbital fog) sector. 
Changes in this external field drives changes in orbital 
occupation on the one hand, and results in a lowering 
of the free energy (of the metallic solution) as a func- 
tion of A a p , stabilizing the metallic solution at a critical 
value of A a/ 3, as found in the calculations. Comparison 
with data for the insulating state requires incorporation 
of the c-axis dimcrization [7] into our DMFT procedure, 
and requires extension of the present approach. Hence, 
we focus only on the metallic state in what follows. 

We now compare our theoretical results to experiment, 
showing how excellent quantitative agreement with the 
full one-electron spectrum in the metallic (R) phase is 
achieved in this picture. The calculations were done at 
Ti = 150 K < T M i (I-phase) and T 2 = 370 K > T M i 
(M- phase). We used parameter values identical to those 
used in Ref. [3]: the LDA bandwidth W ~ 2.0 eV, and 
interactions specified by U ~ 5.0 eV, Jh = 1-0 eV and 
U' ~ (U - 2J H ) = 3.0 eV for the I-phase. To ac- 
count for reduction of the effective interaction due to 
dynamical (selfconsistent) screening processes relevant 
upon the transition to the M-phase, we use reduced val- 
ues of U ~ 4.6 eV, U' ~ 2.6 eV in the M-phase, while Jh, 
being almost unscreened, is kept fixed. As described in 
Ref. [3] , the MIT in our picture results from large changes 
in dynamical transfer of spectral weight from high- to low 
energies following small changes in the renormalized an- 
tiferroelectric distortion in VC>2- 

In Fig. 2, we show the theory-experiment comparison 
for the total (photoemission + X-ray absorption) inten- 
sity, Itotai(u) for the metallic (R) phase. Remarkably, 
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FIG. 2. Comparison of theoretical (LDA +DMFT) result 
for the total one-electron spectral function in the metallic 
phase of VO2 to the experimental results taken from [11] (for 
PES) and from [9] (for XAS). The inset shows a blow-up of 
the PES spectrum. The theoretical result is compared with 
the raw data without subtraction of the 0-2p and surface con- 
tribution. Subtraction of the surface contribution as in [11] 
would lead to even better agreement, as explained in the text. 

the agreement between theory and experiment in the re- 
gion —2.4 < uj < 0.8 eV is excellent. We have not tried 
to fit the data for ui < —2.4 eV and ui > 0.8 eV because 
only part of the 0-2p and e g bands of V projected onto 
the t2g states have been included in our LDA+DMFT 
calculation. We emphasize that even better agreement 
could have been obtained if the subtraction of the sur- 
face contributions would have been carried out without 
taking out the 0-2p "background" in Fig. 1. This would 
have had the effect of increasing the weight in the region 
—0.3 < u < eV and decreasing the weight in the region 
— 1.8 < uj < —0.5 eV [11], further improving the fit. 

Several features of the experimental spectra are read- 
ily interpreted in the light of the theoretical calculation. 
First, most of the low-energy spectral weight actually cor- 
responds to (largely incoherent) electronic states close to 
the minimum of the pseudogap in the spectral function. 
There is no indication of even a strongly renormalized 
FL-like resonance. The DOS at Ep = eV is drastically 
reduced compared to the LDA value [14], in excellent 
agreement with the data, and in full accord with the bad 
metallic resistivity in the R phase. The lower Hubbard 
band feature around —1.2 eV is reproduced very well, 
as is the shoulder-like feature at 0.3 eV in XAS. In [3], 
we found excellent agreement with the optical effective 
mass, m* /m,b — 3.0, a value that is also consistent with 
the observed enhancement of the paramagnetic spin sus- 
ceptibility in the R phase. Thus, these results represent a 
quantitatively accurate description of basic spectroscopic 
and thermodynamic results for the metallic state. 



FIG. 3. Orbital-resolved one electron spectral functions for 
the metallic phase of VO2 (main frame) and for the insulating 
phase without spin dimerization (inset). Only the d yz ± zx or- 
bital DOS crosses Ef in the metallic phase; the d xy DOS 
still shows Mott-Hubbard insulating features, showing the 
"two-fluid" character of the MIT in V0 2 . 

Armed with this agreement, we proceed to analyze the 
MIT in VO2 in more microscopic detail. Fig. 3 shows 
the orbital-resolved DOS for the metallic phase. Clearly, 
only the yz ± zx DOS gives rise to the metallic behav- 
ior; the xy orbital DOS is still gapped at a value close 
to that for the (undimerized, but this is not the real low- 
est state in the / phase) paramagnetic insulating state. 
This observation, constituting an explicit manifestation 
of the two-fluid scenario for the Mott transition in VO2 
driven by multi-orbital correlations, is the central result 
of our analysis. If confirmed by experiment, this result 
would imply concrete evidence of a multi-orbital corre- 
lation driven MIT in VO^- We suggest that orbital re- 
solved electron spectroscopies can resolve this general is- 
sue of importance for transition- metal oxides. Based on 
the above analysis, we expect that only the yz ± zx DOS 
will cross Ep as T is increased above Tmi, while the 
xy DOS will still show insulating behavior. A detailed 
analysis of the anisotropic optical spectra across the MI 
transition driven by strong changes in orbital occupation 
can also verify the proposed scenario [15]. 

The above results are also consistent with the large 
(bad-metallic), linear- in-T resistivity observed in the 
high-T metallic phase of V02- Indeed, within our DMFT 
treatment of local (multi-orbital) electronic correlations, 
the resistivity is expected to show a quadratic temper- 
ature dependence for T < T co h, a lattice FL coherence 
scale, but a linear-in-T dependence above T co h up to high 
temperatures T* ~ J ~ At 2 /U ~ 0.1 eV [16]. For large 
U /t, and near the MIT, the resistivity increases way be- 
yond the Mott limit, and increases without saturation. 
In V0 2 , the MIT occurs at T M i = 340 K, which is much 
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higher than an effective FL coherence scale (which might 
have been visible if Tmi could have been low enough). 
Hence, the low-T FL-like resistivity is not observed, while 
the linear-in-T resistivity up to T* — 900 K without any 
hint of saturation can now be identified with the inco- 
herent local (spin-orbital) moment regime of the metal- 
lic phase in a multi-band Hubbard model. Using this 
analogy, we also see that the optical conductivity will be 
almost completely incoherent in the M phase. 

We have not considered the insulating phase with spin 
dimerization in this communication. To do so requires an 
explicit consideration of the dynamical effects associated 
with the c-axis dimerization. This requires a non-trivial 
extension of DMFT, and we leave it for the future. 

In conclusion, we have shown how the two-fluid pic- 
ture of the Mott transition in VO2 based on the 
LDA+DMFT(IPT) yields excellent quantitative agree- 
ment with the full one-electron spectrum in the metal- 
lic phase. Several features of the bad-metallic trans- 
port arc also understandable in this picture. Orbital 
selective one-electron- and optical spectroscopy can con- 
firm/refute this outstanding issue that, in the final anal- 
ysis, is a manifestation of Mott transitions driven by 
strong, multi-orbital electronic correlations. 
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